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The semileptonic B ■kI'^vi decay is studied starting from a simple quark model which includes 
the influence of the B* pole. To extend the predictions of a nonrelativistic constituent quark model 
from its region of applicability near g^ax = {rnB ~ m^)^ to all values accessible in the physical 
decay, we use a novel multiply-subtracted Omnes dispersion relation, which considerably diminishes 
the form factor dependence on the elastic nB —> nB scattering amplitudes at high energies. By 
comparison to the experimental branching fraction we extract \Vub\ = 0.0034 ± 0.0003 (exp) ± 
0.0007 (theory). To further test our framework, we also study D ~* n and D ^ K decays and find 

O excellent results = 0.80 ± 0.03, 5f?o°^^~''t'''^'^ = 0.079 ± 0.008. In particular for the D tt 

case, we reproduce, with high accuracy, the three-flavor lattice QCD results recently obtained by 
the Fermilab-MILC-HPQCD Collaboration. While for the D —> K case, we successfully describe 
the data for /^(fj'^)//^(0) recently measured by the FOCUS Collaboration. 
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I. INTRODUCTION 

> 

OO . Exclusive semileptonic decays of B-mesons are of great interest, since they can be used to determine the Cabibbo- 
■ Kobayashi-Maskawa (CKM) matrix elements \Vub\ and \Vcb\- In the latter case, heavy quark symmetry greatly 
simplifies the theoretical understanding of the hadronic transition matrix elements and thus the overall theoretical 
uncertainty on the decay process is under control Ij. The measurement of the exclusive semileptonic decay B —^ ttI^vi 
by the CLEO CoUaboration HQ can be used to determine the CKM parameter \ Vub\- This exclusive method provides 
. an important alternative to the extraction of \Vub\ from inclusive measurements of -B ^ Xy^l'^vi. For semileptonic 
decays of charmed or bottom mesons into light mesons there are no flavor symmetries to constrain the hadronic matrix 
— . elements, and as a result, the errors on |T4b| are currently dominated^ by theoretical uncertainties 0. An accurate 
I ■ determination of \Vub\ with well-understood uncertainties remains one of the fundamental priorities for heavy flavor 
Qh, physics. 

^ The transition amplitude for the exclusive semileptonic b —> u decays factorizes into leptonic and hadronic parts. The 
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hadronic matrix elements contain the non-perturbative, strong-interaction effects and have been extensively evaluated 



Deen extensively i 

within different approaches. Thus, several lattice QCD (first in the quenched approximation, B H, la B, IM B i and more 
recently using dynamical configurations 10, ll]), light-cone sum rule (LCSR) [Ullllllli?,"]^ and constituent quark 
model 17, 18, 19^ 20^ 2L> 22ii 2^ 2^ calculations have been carried out in recent years. Each of the above methods 
■ has only a limited range of applicability, namely: LCSR are suitable for describing the low squared momentum transfer 
(g^) region of the form factors, while Lattice QCD, because of the limitation on the magnitude of spatial momentum 
components, provides results only for the high region. Constituent quark models may give the form factors in the 
full range, but they are not closely related to the QCD Lagrangian^ and therefore have input parameters which 
are not directly measurable and might not be of fundamental significance. Thus, it is evident that a combination of 
various methods is required. 

Watson's theorem for the B ttI^vi process allows one to write a dispersion relation for each of the form factors 
entering in the hadronic matrix element. This procedure leads to the so-called Omnes representation |2^, which can 
be used to constrain the q^ dependence of the form factors from the elastic ttB — > iiB scattering amplitudes [2^ . 
In Ref. p^ . once-subtracted dispersion relations were used, and though promising results were found, they suffered 
from sizeable uncertainties because of imprecise knowledge of the ttB ttB phase shifts far from threshold. A recent 
re-analysis of the Omnes representation in this context |28| , has shown that the use of multiply-subtracted dispersion 
relations considerably diminishes the form factor dependence on the elastic ttB — > ttB scattering amplitudes at high 



^ The current best value for | V^^i \ comes from neutrino production of charm off valence d quarks (with the cross section from perturbative 

QCD), rather than from semileptonic D decays. 
^ A rigorous derivation of this approach as an effective theory of QCD in the non-perturbative regime has not been obtained. 
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energies, and more importantly points out that the Omnes representation of the form factors can be used to combine 
predictions from various methods in different regions. 

In this paper we study the semileptonic B ixl^vi decay. We take advantage of the findings of Ref. and 
use the predictions of LCSR calculations &t = Q to extend the predictions of a simple nonrelativistic constituent 
quark model (NRCQM) from its region of applicability (near q^^^ — {ms — w^)^) to all values accessible in the 
physical decay. We also use the available lattice QCD data to test our approach. We use a Monte Carlo procedure to 
find theoretical error bands for the form factors and the decay width. From our estimate of the decay width and the 
branching ratio measurement of Ref. 3] we obtain 

iKblthis work = 0.0034 ± 0.0003 (exp) ±0.0007 (theory) (1) 

To further test this simple framework, we also study the Z? — > tt and D ^ K decays, for which there exist precise 
experimental data and for which the relevant CKM matrix elements (|Vcd| and \Vcs\) are also well known. We find 

/+(0) = 0.63±0.02, /+(0) = 0.79±0.01, = 0.80 ±0.03, -^-^ \ = 0.079 ±0.008 (2) 

/^(O) 0[D" K e + Ve) this work 

The plan of this paper is as follows. After this introduction, we study the semileptonic _B — > tt decay in Sect. ^] 
First we set up the form factor decomposition fSubsect. lII A|) . discuss the valence quark approximation (Subsect. lII B|) 
and the role played by the B* resonance (Subsect. Ill C|) . The Omnes dispersion relation and its application to this 
decay is addressed in Subsect. lII DI and in the Appendix. Finally, in Subsect. Ill El we use our framework to determine 
\Vub\^ paying special attention to estimating the uncertainties of the determination. In Sect. IIIll we study the Z? — > tt 
and D ^ K semileptonic decays and finally in Sect. Hvl we present our conclusions. 

II. SEMILEPTONIC B ^ n DECAYS 

A. Differential Decay Width and Form Factor Decomposition 

Using Lorentz, parity, and time-reversal invariance, the matrix element for the semileptonic — s- ir^l^i/i decay 
can be parametrized in terms of two invariant and dimensionless form factors as'^ 

{7:ip^)\V^\Bips)) = i^PB +P.~ g "'^^', /+(g^) + g^ ""^ ~ /°(g^) (3) 

where q^ = {pB — P-k)^ is the four momentum transfer and ms = 5279.4 MeV and r7i,r = 139.57 MeV are the B^ 
and 7r~ masses, respectively. The physical meaning of the form factors is clear in the helicity basis, in which /+ (/") 
corresponds to a transition amplitude with 1^ (0^) spin-parity quantum numbers in the center of mass of the lepton 
pair. For massless leptons {I — e or the total decay rate is given by 

with g^^^ = {tub - m^f , Gp = 1.16637 x 10"^ GeV"^ and X{q'^) = {m% + ml - q^f - 4m%ml = 4m%\p^\'^, with 
Ptt the pion three-momentum in the B rest frame. 

Measurements of the B'^ lifetime, r^o = (1-536 ± 0.014) x 10^^^ s and of the B'^ — > tt^I^i^i branching fraction, 
Scxp(SO ^ TT-l+vi) = (1.33 ± 0.22) X lO-'' Q lead to 

Tcxp {B° n-l+i^i) = (8.7 ± 1.5) X 10^ s"^ = (5.7 ± 1.0) x 10"" MeV, ^ = e or ^ (5) 

B. Nonrelativistic Constituent Quark Model: Valence Quark Contribution 

Within the spectator approximation, considering only the valence quark contribution and assuming that the B and tt 
mesons are S'-wave quark- antiquark bound states, a NRCQM (with constituent quark masses mb and mi — m„ — md) 



Note that the axial current docs not contribute to transitions between pseudoscalar mesons. 
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predicts |29j: 
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with E^ = y/ml + q"^, Eb^u{k) 



+ k^, nisp the spectator quark mass {md in this case) and (t>,.^\ {k) the 



Fourier transforms of the radial coordinate space B, n meson wave functions, which describe the relative dynamics of 
the quark- antiquark pair"*. 

To evaluate the coordinate space wave function we have used several nonrelativistic quark- antiquark interactions. 
Their general structure is as follows [s^, |^ 
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with <T the spin Pauli matrices, rm the constituent quark masses and 



Xo{mi, rrij 



A 



2mimj 



(7) 



(8) 



The potentials considered differ in the form factors used for the hyperfine terms, the power of the confining term {p — I, 
as suggested by lattice QCD calculations [s^l, or p ~ 2/3 which gives the correct asymptotic Regge trajectories for 
mesons js^l), or the use of a form factor in the one gluon exchange Coulomb potential. All interactions have been 
adjusted to reproduce the fight (tt, p, K, K* , etc.) and heavy-light (D, D* , B, B* , etc.) meson spectra and lead to 
precise predictions for the charmed and bottom baryon (Ac,6, Sc,f), 5]* Sc,&, S* b, ^c,b and il* j,) masses |3lLl33| 



and for the semileptonic A? A+Z and ^? — > ii+Z 



decays. 



Typical NRCQM valence quark predictions for the and /° form factors are depicted in Fig. ^ The ALl 
potential from Ref. 31] has been used^ and for comparison quenched lattice results are also plotted. Preliminary 
unqucnched lattice calculations have been presented recently ^idi , but no significant difference between quenched 
and unquenched calculations is observed [3g, within relatively large statistical errors. In addition, LCSR provide 
accurate and theoretically well founded results in the q^ — region. Thus, we have a LCSR value [13 



LCSR: /+(0) = 0.28 ±0.05 



(9) 



which is also plotted in Fig. ^ 

Fig. n clearly shows the deficiencies of the NRCQM valence quark description of the B ttI'^vi semileptonic 
decay. It fails over the full range of q^ values. Close to Qmax: where the nonrelativistic Schrodinger equation should 
work best, the influence of the B* resonance is clearly visible At the opposite end, close to q^ = 0, where 

\q\ « 2.5 GeV, predictions from a nonrelativistic scheme are clearly not trustworthy. As a result, a value for the width 



^ NRCQM \^ 



l + Vi) = 2.4 



.0032 



X 10 MeV is obtained, which is around a factor of two smaller than the 



CLEO measurement quoted in Eq. jSJ. 



Nonrelativistic Constituent Quark Model: B* Resonance Contribution 



A NRCQM description of the decay process should be feasible in the neighborhood of g^ax- Indeed, this is the case 
for the semileptonic B Dlvi and B D*Wi decays, recently studied in Ref. ,2^ with the same NRCQM as here. 



4 They are normalized to /^^ dkk'^\(j)f^^^ {k)\^ = 1 

^ The sensitivity of the results to the quark-antiquark nonrelativistic interaction will be discussed in detail later. 
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FIG. 1: NRCQM valence quark and form factors from the ALl inter-quark interaction IsH). Lattice data points are 
taken from Refs. (UKQCD) and (APE) and the LCSR estimate at = is given in Ref. (uj. 

The difference here is that, as first pointed out in Ref. ,19], in the chiral limit and as mi, oo, the decay n^l^vi 
should be dominated near zero pion recoil by the effects of the B* resonance, which is quite close to (Z^ax- the 
picture of Ref. [l9j |. the one we will adopt here, the B* contribution plays a role only near (7^a,x: since it is strongly 
suppressed by a soft hadronic vertex. This is in sharp contrast to phenomenological parameterizations of which 
assume it dominates over the full range accessible in the physical decay The B* effects of the type considered 
here are not dual to the valence quark model form factors and must be added as a distinct coherent contribution to 
heavy quark decay near gj^^x [13 ■ ^il^ focus on the /+ form factor, which determines the decay width for massless 
leptons, and we evaluate the contribution to it from the diagram depicted in Fig. |21 It leads to a hadronic amplitude 
(normalizations as in Ref. [2^) 
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with ruB* = 5325 MeV, and /s* and gB'B-K the B* decay constant and the strong B*Btt dimensionless coupling 
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FIG. 2; B* resonance contribution to the form factor for the semileptonic B — > tt decay. 



constant for a virtual B* meson, respectively. On the B* mass shell, the hadron matrix elements /s* (q^ 
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and gB* Btt {q^ = ) = Qb-'Btt reduce to the ordinary B* decay constant and coupling of a pion to B and B* mesons. 
The latter is related, in the heavy quark limit, to g, the coupling of the vector and pseudoscalar heavy-light mesons 
to the pion |37ll3i|^ 



9s^B.= [ ''^Y^' ]il + Oil/m,)) (11) 



From Eq. (|10|) we get 



/polcW ) - ^gB*B.{q ) ^2 _ 2 (12) 



There is no direct experimental determination of gs-BTr, because there is no phase space for the decay B* — > B-k. 
The available experimental results for D* — > Dtt [i|| can be related to gB*BTT, through heavy quark symmetry. There 
is no direct measurements of /b' either. In Ref. [23| we computed, within the same NRCQM approach as the one 
outlined here, both gB-'B-w and and we found a value of 9.1 ± 0.9 GeV for the product of both quantities, which 
appears in f^^]^^ at = m^* . Lattice QCD simulations have measured /b* [39| and gs* b-k |40j | to be 

/b. = 190±30MeV gs-s, = 47 ± 5 ± 8 ^ [5B*s./i3*]Latt-QCD = 8-9 ± 2-2 GeV (13) 

where we have added errors in quadrature. Thus the lattice prediction for the product gB'B-nfB* is in remarkable 
agreement, within 3%, with our NRCQM estimate in ,22,]. In what follows we will use the value and error for the 
product estimated from the lattice data and use the NRCQM of Ref. 0| to determine the dependence of gs* S7r('Z^) 
and fB'iq"^), as we will discuss below. There are other recent estimates for gB'B-n ('4T,'45|) and /s* (^3|), but given 
the existing uncertainties, all of them are compatible with the lattice values quoted in Eq. H13|l . 

As mentioned above, we use the NRCQM framework to estimate the dependence of the product of 
fB*{q^)9B-'BTr{q^)- Since the NRCQM always uses on-shell meson wave functions, all dependence will arise from 
the kinematical factors relating the quark model matrix elements and the hadron form factors. For instance, from 
Eq. (13) of Ref. we find a rather mild q^ dependence 



fB'{q')</q^^fB'V^ (14) 

In the same manner, we use Eqs. (50) and (51) of Ref. to determine the q^ dependence of gB'BTT, setting the 
B-meson four momentum P'^ = {mB,P' ~ 0) and the i?*-meson four momentum = q^ = {niB ~ Ett,P = ^q), 



and off-shell mass given by \/q^, \q\ — y/E^^^-rn^, and determined from q^ as usual (i?^ ~ (m^ +771^ — 5^ )/2mB). 
Thus, finally we evaluate 

/poJ'Z ) = -^GB.{q )^ (15) 

where Gb' (q^) — gB* B-niq^) / dB^Bir is a dimensionless hadronic factor normalized to one at q^ — m^, , which accounts 
for the q^ dependence of B ^ B*tt amplitude. In Fig. 13 we show the influence of the B* resonance within our 
NRCQM and compare our results to those obtained by Isgur and Wise from the gaussian constituent quark model 
of Refs. Our model for the B* contribution compares well to that of Ref. 0, though the latter decreases 

faster owing to the use of a harmonic oscillator basis. The inclusion of f^^i^ clearly improves the simple valence quark 

contribution and leads to a reasonable description of the lattice data from q'^^^ down to q^ values around 15 GeV^. 
The low q^ region is still poorly described within the current model since relativistic corrections there should be large. 

The hadronic amplitude of Eq. H10|l also leads to a small contribution to the form factor. Though it also improves 
the description for the highest q^ values, it is not large enough and it is necessary to consider the influence of the 
lightest 0+ S-resonances (for instance a resonance around 5660 MeV [ilp. 

D. Omnes Representation 

Here we use the Omnes representation of the /+ form factor to combine the NRCQM predictions at high q^ values, 
say above 18 GeV^, with the LCSR result at q^ = 0. In this way we obtain the full q^ dependence of the form factor 



We use the normalization /,r ~ 131 MeV. 
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FIG. 3: Left: The solid line denotes the ALl NRCQM B* pole contribution to (Eq. UlSll 'l. while the dashed line stands 
for the B* contribution to obtained within the gaussian constituent quark model of Refs. [18, 19]. We also plot the G_g* 
hadron factor introduced in Eq. 1151 and the valence quark contribution to depicted in Fig. Right: Valence quark and 
valence quark plus B*— pole (denoted as NRCQM) contributions to We also plot lattice QCD and LCSR /"'' data from 
the same references as in Fig. Q 



and thus can determine the \Vub\ CKM matrix element from the integrated semileptonic width. As shown in the 
Appendix, the (n + l)-subtracted Omnes representation for /+ reads: 



1 r+°° ds 
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(16) 



with < sth = (ws + 'TT-Tr)^ and q^, - ■ ■ , qf^ G ]— oo, sth[- This representation requires as an input the elastic ttB — s- ttB 
phase shift 6{s) in the = 1^ and isospin 1 = 1/2 channel plus the form-factor at (n + 1) g^-values (gg, qf,---, q^) 
below the ttB threshold. 

We would like to stress that from a theoretical point of view the Omnes representation is derived from first principles: 
the well-established Mandelstam hypothesis ji^ of maximum analyticity and Watson's theorem 



f+{s + ie) ^ T{s + ie) 
/+(s-ie) r(.s-ie) 
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Omnes ideas have been used successfully to account for final state interactions in kaon decays ^^'^ -R-^f- 113 j 
a once-subtracted Omnes representation (subtraction point q^ = 0) was applied to the study of semi-leptonic B —> n 



There, however, multiple derivatives evaluated at a single point are used as input instead of subtractions for different values. 
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decays. In the latter work phase shifts were evaluated by solving the Bethe-Salpeter equation in the so-called on- 
shell scheme [49|. with a kernel determined by the direct tree level amplitude from the lowest order heavy meson 
chiral perturbation theory lagrangian : together with the tree diagrams for B* exchange which involve the leading 
interaction with coupling g. Such a model accommodates the B* as a -kB bound state and should acceptably describe 
phase shifts close to threshold. It led to promising results for /+ |23|, but theoretical uncertainties on the form factor 
were not negligible, since to compute the Omnes factor Is (Eq. IjlBfl 'l requires elastic phase-shifts far from threshold®. 
To include the effects of higher resonances on 5{s) requires input of the masses and couplings of such resonances. We 
therefore make many subtractions in the Omnes dispersion relation to suppress the impact of 5{s) at large s. This 
will leave a systematic effect in our results, but this should be less than that coming from the NRCQM plus B* pole 
used as our main input. 

As the number of subtractions increases the integration region relevant in Eq. IjKil) gets reduced and, if this number 
is large enough, only the phase shifts at or near threshold will be needed. Note that close to threshold the p-wave 
phase shift behaves as 

^(s) = nfeTT - p^a H (18) 

where is the number of bound states in the channel (Levinson's theorem H^), p is the ttB center of mass momentum 
and a the corresponding scattering volume. In our case rif, = 1 if we consider the B* as a ttB bound state. Here, we 
will perform a large number of subtractions so that approximating 5{s) ~ t: in Eq. (|16|l will be justified. The Omnes 
factor Is can then be evaluated analytically and we find for < sth 

^ 3=0 

Next we use the above formula to combine the LCSR result at = and those obtained from our NRCQM in the 
high region and presented in the previous section. Thus we have used the NRCQM (valence + pole) predictions 
for five q'^ values ranging from qf^^^ down^ to about 18GeV^: 

{(23.574,4.1373), 
(21.804,2.5821), 
(21.116,2.1969), (20) 
(20.173,1.7916), 
(18.290,1.2591) 

together with the LCSR result of Eq. at = 0. When one uses a large number of subtractions, as is the case here, 
the aj exponents become large and there are huge cancellations (note the normalization condition given in Eq. (|A7|) ) . 
This is the reason why above, and to ensure numerical stability, we have quoted five significant digits for the NRCQM 
input. We are aware that uncertainties are larger than a precision of five digits, and we will carefully take this fact 
into account below. Results are shown in Fig. ^ As can be seen there, we obtain a simultaneous description of both 
lattice data in the high q^ region and the LCSR prediction at = 0. In this way, starting from a nonrelativistic 
valence quark picture of the semileptonic process (Subsect. Ill Bp with all its obvious limitations, we have ended up 
with a realistic description of the relevant form factor for all q^ values accessible in the physical decay. 

A final remark concerns the use of the simplified Omnes representation of Eq. (|19|l instead of the exact one of 
Eq. For instance, if we use five subtractions (we drop the NRCQM point at q"^ — 21.1 GeV^) and the full Omnes 

representation^^ of Eq. H16() . we find tiny differences from the results shown in Fig.^ These differences are negligible 
(below 1%) above lOGeV^, and though larger, still quite small (around 5-7% at most in the 5GeV^ region) below 
lOGeV^ 



* Higher resonance effects on phase shifts cannot be neglected far from threshold. In particular the LCSR result at = hi nts that at 
least an extra = 1~ resonance, located around 6GeV, has to be included in the once-subtracted Omnes relation scheme |27I . 

® Prom Eq. jSJ, we see that the arguments of the meson wave function are \r\ and \l'+ q/2\. For = 18GeV^, half of the transferred 

^ 1 

momentum, \q\/2, is about 0.4—0.5 GeV, which is of the same order as {I '^)^b tt ■ Since the non-relativistic quark-antiquark interactions, 

y(r), have been adjusted to reproduce the meson binding energies, they effectively incorporate some relativistic corrections and hence 

one might expect this effective nonrelativistic framework to provide reasonable meson wave functions for momenta of order ^■ 

This could explain why the NRCQM describes the lattice data (right panel of Fig. |3J from high values of down to values of q^ even 
smaller than 18 GeV Nevertheless, we find it surprising that the nonrelativistic constituent quark model works as well as it does JOj. 
^ We use the model of Ref. to obtain the phase-shifts. 
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FIG. 4: Solid line: /"*" form factor from a 6-subtracted Omnes representation, Eq. 1191 1. Triangles denote the input subtraction 
points of the Omnes dispersion relation (five points obtained from our ALl NRCQM plus the Ref. [H LCSR resuh at = 0). 
Lattice data are from Refs. 0](UKQCD) and (^(APE), and the ALl valence quark plus i?*-pole contribution to /'*' is also 
shown (denoted /nrcqm)- Finally, ±cr lines show the theoretical uncertainty bands on the Omnes form factor inherited from 
the errors in Eq. 1131 and from the quark-antiquark interaction model dependence (see Subsect 111 El for details). 

E. Determination of |Kib|: Error Analysis 

The CKM element \Vub\ can be determined by comparing the experimental decay width (Eq. jSJl) with the result 
of performing the phase space integration of Eq. Q using the form factor /+ determined in the previous subsection. 
Here, we will pay special attention to estimating the theoretical uncertainties. We have two main sources of theoretical 
errors: 

1. Uncertainties in the constituent quark-antiquark nonrelativistic interaction: To estimate those, we will evaluate 
the spread of integrated widths obtained when five different potentials (ALl, AL2, API, AP2 and BD, in the 
notation of Ref. |3l|) are considered. The forms and main characteristics of those potentials were discussed in 
Eqs. 10 and ||HJ). As mentioned in Subsect. Ill Bl all interactions have been adjusted to reproduce the light and 
heavy-light meson spectra and lead to precise predictions for the charmed and bottom baryon masses [sj and 
for the semileptonic A° K'^l~vi and Efl~Di [s^ decays. 

2. Uncertainties on [gB'BjrfB'] cind on the input to the multiply-subtracted Omnes representation: Errors on 
[ffs* Btt/b*] , quoted in Eq. (|13l) . affect the B* pole contribution to /+ (see Eq. (|15|l ') and also induce uncertainties 
in the NRCQM prediction for the five points used as input to the Omnes representation in Eq. (|19|l . Quark- 
antiquark potential uncertainties, discussed in the previous item, also induce uncertainties in the Omnes input. 
The errors on the = data point (LCSR), quoted in Eq. @, should also be taken into account. 

To take these uncertainties into account, we proceed in two steps: 

1. We fix the quark-antiquark potential to the ALl interaction as in all previous subsections. By means of a Monte 
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Carlo simulation, we generate a total of 1000 ([ffB*B7i-/B*]Latt-QCD i /'''(O)lcsr) pairs^^ from an uncorrelated 
two dimensional gaussian distribution, with central values and standard deviations taken from Eqs. H13|l and iPJ, 
respectively. For each of the 1000 pairs we build up the six points that we use in our Omnes scheme (Eq. (|19|) ) 
and thus find 1000 different determinations of /"*" over the whole range accessible in the B ^ tt decay. For 
each value of q^, we discard the highest and lowest 16% of the values obtained for the form factor, to leave a 
68% confidence level band which forms part of the theoretical uncertainty shown in Fig. ^ Since the output 
distributions are not gaussian in general, this accounts for possible skewness. 

Performing the phase space integration for each of the 1000 form factor samples and again discarding the highest 
and lowest 16% of the values, we find 

^ ^ - (0.50^°;}^) X 10-^MeV (21) 



\Vul 



12 ~ v"->j"-o.io; 



2. We fix the ([5B*S7r/B*]Latt-QCD i/^(0)lcsr) pair to their central values and compute the decay width with 
each of the five quark- antiquark interactions discussed above. From the spread of output values, we find 

— ■ ^ = (0.50 ± 0.15) X 10-* MeV (22) 

\Vubr 

Adding both sources of error in quadrature, we get 

T{B" TT-l+vi) 

and comparing to the measurement of the width in Eq. (0 we find 

|K;,|this work = 0.0034 ± 0.0003 (exp) ± 0.0007 (theory) (24) 

The CLEO Collaboration Q obtains from studies of the B*^ tt^I^vi branching fraction and q^ distributions, using 
LCSR for < g2 < 16 GeV^ and lattice QCD for 16 GcV^ <q^ < g^ax, 

+0.0006 



(0.50 ± 0.20) X 10"» McV (23) 



iKblcLEO = 0.0032 ± 0.0003 (exp)+^:^^^^ (theory) (25) 

We see that both determinations of \Vub\ are in an excellent agreement and that in both cases the error is dominated 
by uncertainties in the theoretical treatment. We have also calculated partially integrated branching ratios. 



S*°*p^1(B0 ^ TT~l+vi) fi2 ^ 2 dr 

" d? 



Biqi <q'< qi) = dq^ (26) 



Theoretical uncertainties partially cancel in the ratio j^i dq^-§^/^- Our results are compiled in the TableQ] There it 
can be seen that they compare reasonably well with those quoted in Ref. 3] . 



CLEO \3] 
This work 



e(0 < < 8 GeV2)/10-^ B(8 < < 16 GeV^)/10-'^ B{q^ > 16 GeV2)/10 



4.3 ±1.2 6.5 ±1.3 2.5 ±1.0 

4.3 ± 0.7 [exp] ± 1.2 [theory] 4.1 ± 0.7 [exp] ± 0.4 [theory] 4.9 ± 0.8 [exp] ± 1.2 [theory] 



TABLE I: Partially integrated branching ratios (see Eq. 12611 '). 

Finally, at each value of q^ we also compute the spread of values obtained for the /+ form factor when the five 
different quark-antiquark interactions are used. This procedure gives us a further theoretical error on /^(g^) at fixed q^ 
and by adding it in quadrature to that obtained previously from uncertainties on the {[gs' STr/s'lLatt-QCD ' /^(O)lcsr) 
pair, we determine the theoretical error bands shown in Fig. 0] 



We have checked that the errors quoted in the following are already stable when 500 event simulations are performed. 
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III. SEMILEPTONIC D ^ tt AND K DECAYS 

As a further test of our predictions for the B —t t: semileptonic process, we present results for the D tt and 
D K decays for which there are precise experimental data [j]: 

Bc^piD" TT-e+j^e) = (3.6 ± 0.6) x lO^^, T^^p{D^ i^'e+v^) = (5.8 ± 1.0) x IQ-^^ MeV (27) 
B^^p{D^ R-e+Ve) = (3.58 ± 0.18) X 10-2, rcxp(£'° -> i^^e+i/e) = (57 ± 3) x lO^^^ MeV, (28) 

with life time r^o = (410.3 ± 1.5) x lO^^^ g ^nd \Vcd\ = 0.224 ± 0.003, \Vcs\ ^ 0.9737 ± 0.0007. 

In the last two years there has been renewed interest in these decays. The first three-flavor lattice QCD results [S^ 
have appeared, superseding the old quenched ones [111113, and the BES |^] and CLEO |E3 collaborations have 
new measurements of the branching ratios 

BES : S(D" ^ TT^e+i^e) = (3.3 ± 1.3) x 10-^ ^(1?" ^ if-e+i/e) = (3.8 ± 0.5) x 10"^ (29) 

Both collaborations have also determined the form factor at = 

BES : /+(0) = 0.73 ±0.14 ±0.06, /^(O) = 0.78 ± 0.04 ± 0.03, ^^^^ = 0.93 ± 0.19 ± 0.07 (31) 

CLEO:^^^ ^ 0.86 ±0.07t[5;[5^± 0.01 (32) 

In the following wc will apply the NRCQM developed for the B ^ tt decay to the description of these Z3-mcson 
semileptonic transitions. All formulae of Sect. can be used here with the obvious replacements: B D, B* ^ D* 
for the D° tt~1'^vi process, and B ^ D, n ^ K, B* ^ D*^ for the Z?° K~l+vi process. We will use 
■niDO = 1864.6 MeV, mo- = 2010 MeV, rrio'^ = 2112.1 MeV and rriK- = 493.68 MeV. 

A. D -kIvi 

Since there is phase space for the D* — > Z?7r decay to occur, the qu'Dt! hadronic constant has been experimentally 
measured (CLEO 54]) 

g^,.^^ = 17.9 ±0.3 ± 1.9 (33) 
Taking fo- = (234 ± 20) MeV from Ref. HI, we find^^ 

[ffi5.i5./i5.]Exp-Latt = 4.2 ± 0.6 GeV (34) 

where we have added errors in quadrature. In Ref. [2l| and using the same set of NRCQM's, we found a value of 
4.9 ± 0.5 GeV for the above product, in reasonable agreement with Eq. (|34|l . The value quoted in Eq. H34I) determines 
the £'*-pole contribution, above = 0, to /"*" and adding it to the valence quark contribution we obtain the results 
shown in Fig. [3 We find excellent agreement between our description of the form factor and that provided by the 
unquenched lattice simulation of Ref. [Hl|. As can be seen in the figure, the Z)*-pole contribution is dominant above 
= 1.5 GeV^ and it remains sizeable down to 0.5 GeV^. We do not see the need to Omnes improve the NRCQM 
description of the decay since it is quite good for the whole q^ range. On the other hand, we see that the pion energy 
ranges from m^r up to about 1 GeV, which was also the maximum value for E^^ in the five NRCQM data-points used in 
the subtracted Omnes representation for the semileptonic B — > tt decay depicted in Fig. 0] This reinforces our belief 
in the reliability of our determination of \Vub\ presented in Subsect. Ill El Considering our theoretical uncertainties 
(errors on Eq. H34|) and the spread of results obtained when different quark-antiquark interactions are considered) 
together with the experimental uncertainties on |Vcrf| quoted above, we find 

Tthis work(i?° ^ TT-e+i/e) = ^5.2 ± 0.1 (cxp :|T4<j|) ± 0.5 (theory)! x IQ-^^ MeV (35) 



Note that the lattice QCD simulation of Ref. ^5E\ measured go'Dn = 18.8 i 2.3 ± 2.0 in good agreement with Eq. 1331 . 
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FIG. 5: Valence quark (dotted line) and valence quark plus Z)*-pole (solid line denoted /nrcqm) contributions to for 
D — > TT semileptonic decay. In both cases the ALl quark-antiquark interaction has been used. Triangles (Q) and circles 
(113) stand for the lattice QCD quenched results obtained by the APE and UKQCD Collaborations, respectively. We also 
plot the three-flavor lattice QCD results [13| from the Fermilab-MILC-HPQCD Collaboration (FMH, diamonds), the best fit 
(dash-dotted line) to this latter set of data-points (Eq. (5) of Ref. [Holl'). and the determination of /+ at g'^ = (square) by the 
BES Collaboration 52] . Finally, the ±ct lines stand for the theoretical uncertainty bands, inherited from the errors in Eq. 13411 
and from the quark-antiquark interaction model dependence. 



in good agreement with Eq. 1)2 7|l . We also obtain 



/+(0) =0.63 ±0.02 



(36) 



compatible within errors with both the BES (0.73 ± 0.15) and the Fermilab-MILC-HPQCD (0.64 ± 0.07) results. 

Finally in the left plot of Fig.EJ we compare the NRCQM predictions for the ratio f^{q^)/ f^{0) with a pole form 
recently fitted to data by the FOCUS Collaboration j56l |. 



B. KWi 

Since there is no phase space for the D* — > DK decay, we will estimate the go'DK coupling from the value quoted 
for go-'DTT in Eq- H33() . The parameter g defined in Eq. 111(1 describes the strong coupling of charmed mesons as well 
as of beauty mesons to the members of the octet of light pseudoscalars. We will assume fiavor SU(3) symmetry for 
this basic quantity in the heavy quark chiral effective theory, and thus we will use [H^ 

gD-DK~ 7 ~gD'Dvr , — 15.3 ± 1.6 (37) 

JK ^/mo* JK 

where we have taken Jk/ fir ~ 1.2 from Ref. [H^, and have kept some SU(3) flavor breaking terms in the masses of 
the charmed vector mesons and in the kaon decay constant. Taking fw — (254 ± 15) MeV from Ref. [3^, we find 

[5D|^^/^dsu(3)^Latt =3-9±0.5GeV (38) 
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[GeV2] [GeV^] 

FIG. 6: NRCQM predictions for the ratio /+(g^)//+(0) for both D ^ tt (left) and D ^ K (right) semileptonic decays. For 
comparison we also plot experimental results from the FOCUS Collaboration fs^ : pole fit for the D n decay (mpoie = 
I.OI^q 'J^ GeV) and direct measurements of the form factor for different values, in the D ^ K case. In this latter case, 
the ±(T lines stand for our theoretical uncertainty bands, inherited from the errors in Eq. 1341 and from the quark-antiquark 
interaction model dependence. 



Our results for this decay are shown in Fig. Several comments are in order: 

1. We find reasonable agreement with the three-flavor lattice QCD results [S^ up to kaon energies of the order 
of 1 GeV, which cover the whole range accessible in the physical decay. Discrepancies with lattice data are 
now more sizeable than in the _D — s- tt case and lattice unquenched data favor values for [ff-Dj-D-ffZ-Dj ] su(3)_Latt 
smaller than the one used in our calculation (see Eq. 1123). Theoretical errors for /+ arc, in this case, mostly 
due to the uncertainties on [ff-Djr'-R'/-D'] gu(3)_Lj^j^- Nevertheless, we would like to point out that uncertainties 
on the value of gn^DK might be larger that those quoted in Eq. 137|) . since flavor SU(3) corrections to the 
relation cjd'DK ~ ^g^mjjmi:)' / Jk could be large (ms/mc ^ „/mc). 

2. The contribution of the vector resonance is less important than in the B ^ tt and D ^ it decays, since the D* 
is located relatively far from vymax- 

3. Our predictions for /+ at negative values of , which do not enter into the phase space integral, suffer from 
larger uncertainties, since in that region the transferred momentum is larger than 1 GeV and, as for the B t: 
case, relativistic effects could became important. One could Omnes improve the NRCQM to achieve a better 
description of the form factor in the negative q^ region. 

In the right plot of Fig. El we compare the NRCQM predictions for the ratio /^(q^)//^(0) with recently measured 
data from the FOCUS Collaboration ,52J and find satisfactory and reassuring agreement. 
For the integrated width, we find 



this 



i,{D° ^ R-e-^Ve) = 66 ± 3 (theory) x 10"^^ MeV (39) 



which is about two standard deviations higher than the value quoted in Eq. (|28|) . 
Eqs. H27() and (|28|l lead to (adding errors in quadratures) 



mi^o^^-e.^j) - 0-101 ±0-01^ (40) 
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FIG. 7: Valence quark (dotted line) and valence quark plus D*-pole (solid line denoted /nr^cqm) contributions to /"*" for D ^ K 
semileptonic decay. In both cases the ALl quark-antiquark interaction has been used. Circles (j5L]) and diamonds (HOl) stand 
for the lattice QCD quenched and unquenched results obtained by the UKQCD and the Fermilab-MILC-HPQCD Collaboration 
(labelled FMH), respectively. We also plot the best fit (dash-dotted line) to this latter set of data-points (Eq. (5) of Ref. |HOI), 
and the determination of /"*" at g'^ = (squared) by the BES Collaboration [H^l . Finally, the ±a lines stand for the theoretical 
uncertainty bands, inherited from the errors in Eq. 13411 and from the quark-antiquark interaction model dependence. 



which turns out to be a bit higher, though compatible within errors, than the recent CLEO determination quoted in 
Eq. H3Uf) . For this ratio of branching fractions we find 



0.079 ±0.008 



this work 



in excellent agreement with the CLEO measurement. We also find 

/+(0) = 0.79 ±0.01, = 0.80 ± 0.03 



(41) 



(42) 



which compare well to the recent experimental measurements in Eqs. H31() and (|32|l . 



IV. CONCLUDING REMARKS 



We have shown the limitations of a valence quark model to describe the B ^ tt, D ^ n and D ^ K semileptonic 
decays. As a first correction, we have included in each case the heavy-light vector resonance pole contribution. 
For the semileptonic B ^ tt decay, the inclusion of the B* degree of freedom provides a realistic dependence 
of the relevant form factor, /+, from q^^^ down to around 18GeV^. We then use a multiply-subtracted Omnes 
dispersion relation, which considerably diminishes the form factor dependence on the elastic ttB ttB scattering 
amplitudes at high energies, to combine LCSR results at = with NRCQM predictions in the high q^ region. 
As a result we have been able to predict the /+ form factor for all q^ values accessible in the physical decay. We 
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have used a Monte Carlo procedure and analyzed the predictions of five difi^erent quark-antiquark interactions to 
determine theoretical error bands for form factors and the decay width. This has allowed us to extract from the 
measured branching fraction the value \Vub\ = 0.0034 ± 0.0003 (exp) ± 0.0007 (theory) in excellent agreement with 
the CLEO Collaboration determination of Ref. 3J. For the D — > tt semileptonic decay we have found excellent 
agreement between our model calculation (valence quark plus D*-pole contributions) of /"*" and the one obtained by 
the unquenched lattice simulation of Ref. [5fil |. We found no need to Omnes-improve our calculation in this case. 
Our results T{D° ir'e+iye) = [5.2 ± 0.1 (exp :|Kd|) ± 0.5 (theory)] x IQ-^^ MeV and /+(0) = 0.63 ± 0.02 are 
in good agreement with experimental data. Finally for the D K semileptonic decay we find good agreement, 
in the physical region, between our model calculation (valence quark plus D*-pole contributions) of /+ and lattice 
data, from UKQCD jsj and Fermilab-MILC-HPQCD 5d\, and also with recent measurements from the FOCUS 
Collaboration 56]. Again our results B{D^ 7r-e+z/e)/S(L'° -> R-e+v^) = 0.079 ± 0.008, f^{<S) = 0.79 ± 0.01 and 
/^(0)//^(0) — 0.80 ± 0.03 are in good agreement with experimental determinations by the CLEO and BES [s^ 
Collaborations. 
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APPENDIX A: MULTIPLY SUBTRACTED OMNES DISPERSION RELATION 



Let the form factor^'^ /^(s) be analytic on the complex s plane (Mandelstam's hypothesis ji^ of maximum an- 
alyticity) except for a cut L = [sth — {mB + "^7r)^,+oo[ along the real positive s axis, as demanded by Watson's 
theorem ^A^. For real values s < sth the form factor is real which implies that the values of the form factor above 
and below the cut are complex conjugates of each other: f~^{s + ie) — f~^{s — ie)*. For s > Sth, the form factor has a 
discontinuity across the cut and develops an imaginary part /+(s + ie) — /+(s — ie) = 2ilm/(s + ie). Cauchy's theorem 
implies that /^(s) can be written as a dispersive integral along the cut and performing one subtraction at sq < sth 
one gets: 



/+(s)-/+(.so) + 



So 



dx Ini/+(x) 

X — So X — s 



s^L, So < Sth 



(Al) 



Depending on the asymptotic behavior of /^(s) at the extremes of the cut L, more subtractions may be needed to 
make the integral convergent. For the time being, let us assume that one subtraction is sufficient. The well known 
Omnes solution for the above dispersive representation is !2^ : 



0(s) -/+(so)exp 



_so /-+°° dx 8(x) 
■ J,^^ x-sox- s 



S ^ L, Sq < Sth 



(A2) 



with S{s) the elastic ttB ttB phase shift^** in the = 1 and isospin 1/2 channel (see Eq. (^TJ). 0{s) gives the 
physical form- factor since. 



1. For s > Sth, we have 



0(s±ie) = /+(so)exp 



s - So 



dx S{x) 6{s) 

V I ± ITT 

X — Sq X ~ S S — So 



The discussion below can be trivially generalized to any scattering amplitude or form factor with definite total angular momentum and 
isospin quantum numbers. 

Obviously S{s) has to be defined as a continuous function of s. 



/+(so)exp 



So 



V 



dx 5{x) 

sth X-SqX- S 
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(A3) 



where V stands for principal part of the integral. Thus we have 0{s + ie) = 0{s — ie)*, the function O is real 
for s < Sth and has neither poles nor cuts, except for that required by Watson's theorem: L = [sth,+oo[. The 
discontinuity across this cut is given by 0{s + ie) — 0{s~ie) — 2iImO(s + ie) and by construction O{so) — /^(sq). 
Thus, both f~^{s) and 0{s) satisfy the same dispersion relation (Eq. IjAip ) and therefore both functions can 
differ at most by a polynomial with real coefficients, which should vanish at s = sq. But, this polynomial is zero 
since: 



2. The function 0{s) satisfies Watson's theorem: 



(A4) 



Performing n + 1 subtractions, one can can produce a rank n + 2 polynomial in the denominator of the dispersive 
integral of Eq. (|A1() . Indeed, for s ^ L 



(s - sq){s - Si) • • • (s - s„) 



+ 00 



dx 



Im/+(a;) 



TI" Jsth [x - So){x - Si) ■ ■ ■ {x - Sn) X~S 

with the rank n polynomial P,i(s) determined by the n + 1 equations Pn{si) = .f^{si), i = 0, 1, • • • , n 



So, ■ • • , Sn < Sth (A5) 



j), aj{s) = 



n 



S - Sfc 



j#fc=0 



Sk 



(A6) 



Note that aj(s) are rank n polynomials, which satisfy 



On the other hand for s > sth, we have from Eq. (|17|l 



log/+(s + ie) - log/+(s - ie) = log 



/+(s + ie) 
/+(s-ie) 



2i(5(s) ^ 2ilm [log/+(s + ie)] 



(A7) 



(A8) 



Thus in analogy to Eq. (|A5|) . assuming that the form factor does not vanish in C — {sth}^^, or neglecting the contribution 
from the log cut if it has a finite branch point different from sth, we can write 



log/+Gs) =P„(s) + 



(s - so)(s - Sl) • • • (s - s„) 



dx 



5{x) 



{X - Sa){x - Sl)- - - {x - Sn) . 



s(^L (A9) 



with 



P„(s)=^a,(s)log/+(s, 



(AlO) 



From the above equation one readily finds the (n + l)-subtracted Omnes representation given in Eq. H16() . 

Finally, we would like to draw the reader's attention to a subtle point. In Eq. HA1|I we have assumed that /+ has 
no poles. However we know that if the scattering amplitude has a pole at sr = AI^ — iMj^Tn on its second Riemann 
sheet (resonance) or on the physical sheet (bound state with F/; = 0+ and < Sth-), it might show up as a pole 
in the complex plane of /+ (see Eq. (|12|l '). On the other hand, the 5- matrix depends on exp (2i(5), and thus one has 
the freedom to add factors of m7r, for m an integer, to the phase-shift without modifying the S'-matrix. However, the 
Omnes representation of the form factor will definitely depend on the specific value chosen for the integer m. 



Note that, we have already treated s branch point. 
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To fix this ambiguity, we will assume that at threshold the phase shift should be (5(sth) = nbir, where rif, is the 
number of bound states in the channel, while 6{oo) = kn, where k is the number of zeros of the scattering amplitude 
on the physical sheet (this is Levinson's theorem This choice for the phase shifts also takes into account the 

existence of poles in the scattering matrix. We demonstrate with a simple example in which a p-wave T-matrix is 
proportional to (s — Sth)/(s — Af_R, + ^MbTr)- The phase shift is given by 



S{s) — TT + Arctan 



-MrTr 

''-Ml 



s > sth (All) 



with Arctan G [— tt, 7r[. This satisfies d{oo) = n and, if MrTr <C |s — M|,|, it also leads to 5(sth) = tt or for a bound 
state or resonance respectively, in accordance with Levinson's theorem. For simplicity, let us also assume Tr <C Mr. 
In this circumstance we can approximate 

S{s) ^tt[1~ H{Ml - s)] = nH{s - (A12) 

where H{ ) is the step function. Since 

s-sa r+°° dx 5{x) ^ s ) '^^ ^ _ fMax(sth,M2) - So' 



we find that the Omnes solution from a once-subtracted dispersion relation, Eq. (|A2p . reads 

Ois) ^ riso)^pp^^^$^ (A14) 
Max(sth,M|,) - s 

It has a pole at s = for a resonance or at s = Sth in the case of a bound state. For the case of a resonance with 
a finite width, the pole will move to s = + iM^jr^j. On the other hand, for a bound state, going beyond the 
approximation 6{s) « tt (see Eq. H18|l ^. the form factor will be sensitive to the exact position of the pole (s = Af|.), 
since the effective range parameters (scattering volume, . . . ) will depend on Mjj. 

These conclusions can easily be generalized when a multiply-subtracted Omnes dispersion relation is used. 
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